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Abstract 
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Green's functions is presented. Its effectiveness is demonstrated through its application to singular 
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tivistic limits of quantum field theory is highlighted. Not surprisingly, the analogue of ultraviolet 
regularization is required for the interpretation of these singular problems. 
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I. INTRODUCTION 



It has been generally recognized that the treatment of singular problems in quantum 
mechanics is notoriously difficult and plagued with ambiguities [1]. A complete solution 
of a typical quantum-mechanical problem with a singular potential requires the analysis 
of the scattering and bound-state sectors supplemented by a prescription, which can be 
conveniently formulated in the field-theory language of regularization and renormalization. 
Such a prescription establishes the necessary connection between the physics of a particular 
system and the corresponding renormalization of an ill-defined problem. As this paper 
demonstrates, much work is still needed in formulating and interpreting singular problems 
in quantum physics, particularly within the path-integral framework. 

It is also generally accepted that our current understanding of bound states in quan- 
tum field theory — including relativistic effects and radiative corrections — is still unsatisfac- 
tory [2, 3]. Conceptually, the source of this deficiency lies in the apparent contradiction 
between the usual perturbative framework of quantum field theory and the intrinsic nonper- 
turbative nature of bound states, which manifest themselves as poles of the S-matrix and 
can only be generated by summing an infinite set of Feynman diagrams. Therefore, one 
is forced to face simultaneously a plethora of technical difficulties, including issues of regu- 
larization and renormalization, nonrelativistic reductions, model dependence, and stability 
and controllability of the approximation scheme. Following the seminal work of Bethe and 
Salpeter [4] many techniques and ideas have been developed, including numerical schemes, 
and have been applied to both QED and QCD with various degrees of success. More recently, 
the use of effective Lagrangians has advanced the field considerably [5], and extensions to 
other fields of medium- and high-energy physics are underway. 

In this paper we begin a program centered on the use of path integrals for a systematic 
study of singular problems in quantum physics and of the connection between field theory 
and quantum mechanics vis-a-vis the description of bound states. One of the main purposes 
of this approach is to formulate a complementary point of view that may shed light on the 
difficult problems posed by bound states in quantum field theory. It is to be expected that a 
more thorough understanding of these questions would eventually have a substantial impact 
on the development of modern nuclear physics within the chiral Lagrangian approach [5], 
as well as on atomic and molecular physics [1]. Even though the use of path integrals to 
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study bound states is hardly new, it has mainly been implemented within the framework 
of quantum mechanics and, until recently, only for regular potentials [6, 7]. In this paper, 
as we take the first steps in formulating this program, we are confronted with the problem 
of describing singular potentials beyond the Coulomb problem. In particular, we will con- 
centrate upon the path-integral treatment of the very important case of the inverse square 
potential, 1/r 2 , the first-one with integer powers beyond the Coulomb case. 

The quantum-mechanical inverse square potential has been extensively dealt with before 
via the Schrodinger equation [8, 9, 10, 11], and to a lesser extent with path integrals [12, 
13, 14, 15, 16]. A recent proposal for a comprehensive treatment of this problem was 
advanced within the Schrodinger-equation approach, properly combined with regularization 
and renormalization concepts borrowed from quantum field theory. This was done first 
for the one-dimensional case with a real-space regulator [9]; later, it was generalized to a 
complete analysis of the D-dimensional case using different regularization schemes [10, 11]. 
The current relevance of this problem stands out not only from its field-theory treatment but 
also from its recent application to the interaction between an electron and polar molecule, 
which was shown to be a simple manifestation of a quantum anomaly [17]. 

The main goals of the present paper are: 

(i) To introduce a general framework for the computation of bound states and other 
"nonperturbative" results within the path-integral approach to quantum physics (Sec. II). 

(ii) To reformulate the problem of the inverse square potential and reproduce the results 
of Ref. [10] within the path-integral framework (Sees. IV- VI). 

(iii) To generalize the inverse square potential by the inclusion of delta-function interac- 
tions, with miscellaneous applications (Sec. V). 

(iv) To highlight the emergence of singular potentials within the effective-field theory 
program and to sketch possible generalizations in quantum field theory (Sec. VII). 

The treatment of more singular potentials and the development of the full-fledged 
quantum-field-theory case is in progress and will be reported elsewhere. 
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II. GENERAL FRAMEWORK 



A. Quantum-Mechanical Propagator 

Consider now a particle of mass M subject to an interaction potential V(r, t) in D- 
dimensions. Its physics is completely described by the quantum-mechanical propagator or 
transition amplitude 



K D (v",v';t",t') 



Texp 



h 



t" 

Hdt 



(1) 



K D (r",v';t",t') = ^ Vr(t) exp jis[r(f)] (r", r'; t", t') } (2) 



where T is the time-ordering operator and H the Hamiltonian. The propagator admits the 
usual representation 

-r(t")=r" 
'r(t')=r' 

as a path integral, in which S [r(t)] (r", r'; t", t') is the classical action functional associated 
with "paths" r(t) that connect the end points r(i') = r', r(t") = r" . In anticipation of the 
derivations of Sec. Ill, we have carefully defined all the relevant variables in the functional 
dependence of S [r(t)] (r", r'; t", £'); specifically, the action is (i) a functional of r(t) and (ii) 
a function of the end-point variables (r", r'; t", t'), for every chosen r(t). 

Equation (2) is the primitive construct from which we will advance a general technique for 
the evaluation of the physics of singular potentials. Furthermore, we would like to explicitly 
show how to supplement this technique with regularization a la field theory for the derivation 
of the renormalized physics of the inverse square potential — by itself and also combined with 
a delta-function interaction. 

A careful analysis shows that one must exercise proper caution when evaluating Eq. (2) 
in non-Cartesian coordinates. For example, when considering central potentials and trans- 
forming into hyperspherical coordinates in D dimensions, one may need to consider extra 
terms of order h 2 in the action. These extra terms arise whenever nonlinear transformations 
are implemented, a fact that has been established in quantum mechanics [18] and quantum 
field theory [19]. However, these technical difficulties can be altogether circumvented by 
adopting a time-sliced expression in Cartesian coordinates as the starting point in Eq. (2) 
and transforming coordinates before taking the continuum limit. Therefore, we introduce a 
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partition of the time T — t" — t' into N equal discrete time units, 

V (N) {t',t"}: {t = t',t 1 ,...,t N - 1 ,t N = t"} , 

which yields a time lattice of (N + 1) equidistant points 

t" - t' 



tj = t' + je 



N 



(j = 0,...,N + l), 



(3) 



(4) 



with the end points fixed as in Eq. (3). This procedure should be followed by taking the 
limit N — > oo only after all other calculations are completed. Correspondingly, in Eq. (2), 
the action is replaced by its discrete counterpart 



N-l 



3=0 



2e 



(5) 



in which rj = r(tj), while the end points are r = r' and r N = r". It should be noticed 
that the variable N has been isolated as a superscript in order to emphasize its distinct 
nature in the limit S (N) [r 1; . . . r^-i] (r", r'; t", t') -> S [r{t)\ (r", r'; t", If), with the variables 
{i - !, . . .r^v-i} asymptotically generating a "path" r(t). Then, Eq. (2) can be interpreted as 
a shorthand for the time-sliced expression 

rr(t")=r" 



K D (r",r';t",f) = lim 



(*)= 



vWr(t) exp \ j_S^ [n, . . . r^] (r", r'; t", f ) \ , (6) 



where the integration and measure symbols stand for the formal (N — l)-fold integral oper- 
ation 



r r(t")=r" [N-l . 

/ V^v(t) = [C(M,e,h,D)] N JT / d D r k 

Jr(t')=T' Lfc = l^ RD 



(7) 



with 



D/2 



(8) 



As is well known, the propagator Kd(t", r'; t", t') solves the initial- value problem for the 
Schrodinger equation, permitting the computation of the wave function \I/(r, t) at any time 
t, given a particular-time value \l/(r, t ). A related quantity is the retarded (causal) Green's 
function G D (r", r'; t" , t'), which solves the Schrodinger equation with an additional arbitrary 
inhomogeneous term, and which is related to the propagator by 



G D (r", r'; f, t') = 6{t" - f) K D (r", r'; t", f ) , 



(9) 
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where 9(T) stands for the Heaviside function. 

In this paper we will consider time-independent singular potentials, with special emphasis 
on the inverse square potential and the delta-function interaction. These particular cases 
are included in the class of time-independent Hamiltonians, for which the lattice action 

K . . . r N ^] (r", r'; t", f) = S™ [r 1; . . . r^] (r", r'; t" - t>) , (10) 

as well as its continuum limit 

S[v(t)](r",r';t",t') S [r(t)\ (v" ,r';t" - f ) , (11) 

and the propagator 

K D (v",v';t",t') = K D (v",v';t"-t') (12) 

are functions of the end-point times t' and t" only through the difference T = t" — t'. 
Under this assumption, it is also possible to specialize to "paths" connecting the end points 
r(0) = r',r(T) = r". Then, the Fourier transform of Gd(y", r'; T) with respect to T defines 
the corresponding energy Green's function 

G D (v",v';E) = ^ / dTe lET ' h G D {v",v'-T) 



ih 

1 

ih 



4 / dTe iET ' h K D {v\r'-T) (13) 
{ " Jo 

1 r°° rr(T)=r" /■ \ 

= - / dT Vv(t) exp 1 {S [r(t)\ (r", r'; T) + ET}) . (14) 

in Jo Jv(o)=r' \n J 

The energy Green's function is a convenient tool for spectral analyses because it satisfies 
the operator relation 

-l 



G d (t",t';E) = (v" 



(E-H + iej 



(15) 



where ie is a small positive imaginary part — this follows by combining Eqs. (1) and (13). In 
particular, it admits a spectral representation, 

' r '^-Z. E -E n + ie + da E-E a + ie ' (16) 

which identifies the bound states as poles and the scattering states as branch cuts, as well 
as the corresponding stationary wave functions [6]. For example, for a nondegenerate bound 
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state of energy E n , the wave function ip n (r) = |-?/> n (r)|e l<5n ( r ) can be uniquely extracted, up 
to an arbitrary global phase factor, from 

|^ n (r)| = |Res{G(r,r;£ n }| 1/2 (17) 

and 

*„(r) = -<lJ Res{G(r,r ;E n )} 1 

\ [Res {G(r, r; E n )} Res {G(r , r ; E n )}] /2 J 

where r is a convenient reference point and Res{/(z )} represents the residue of the function 

f(z) at the pole z = zq. 

In what follows our arguments and derivations will be based on Eqs. (6) and (14), as well 

as on their counterparts in hyperspherical coordinates, which we will consider next. 



B. Central Interactions 



For the all-important class of central potentials, Eq. (6) can be conveniently rewritten 
in hyperspherical polar coordinates [10, 11, 20]. This can be accomplished by separation 
of variables, a procedure that can be systematically implemented by repeated use of the 
addition formula [21] e izcos ^ = (iz/2)- u T(v) J2Zo( l + v)Ii+ v (iz)C$ v) (cosif>), where I p (x) is 
the modified Bessel function of the first kind and order p, c\ v \x) is a Gegenbauer poly- 
nomial, and v = D jl — lisa related dimensionality parameter. In turn, the Gegenbauer 
polynomials can be resolved into hyperspherical harmonics Yj m (f2), labeled with the IP- 
dimensional angular momentum quantum numbers I and m, where the latter has a multi- 
plicity di = (21 + D - 2){l + D - 3) \/l\(D - 2)! [20]. As a result, the propagator admits the 
partial-wave expansions 

K D {v",v'-T) = ^ (r V)- (D " 1)/2 ^{l + v)Ci v \cos^v)K l+v {r"y-T) 

1=0 

oo di 

= (r V)- (D - 1)/2 E Y lm (QT)Y^(ir) K l+U (r", r>; T) , (19) 

1=0 m=l 

where cosipn"^' = ?" • f' and the radial prefactors are conveniently written so that, for each 
angular momentum channel I, the radial propagator Ki +U (r", r'\ T) satisfies the composition 
property 

poo 

K l+V (r", r'; t" - t>) = / dr K l+V (r", r; t" - t)K l+v (r, r'; t - t') . (20) 

Jo 
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Moreover, the radial propagator in Eq. (19) can be explicitly evaluated by means of the 
path-integral expression 



K l+u (r",r';T) = Jim 



(JL.Y" II / 



N/2 N-l 



dr k 



x exp{j_RW[r 1 ,...r N _ 1 ](r",r';T) 



(21) 



where the radial action is 
In Eq. (21) a radial functional weight 



N-l 



j=0 



M(r J+1 - rj y 
2e 



eV( rj ) 



(22) 



N-l 



A N + l\r 2 } = II [v^ _ *'W*i)] 

j=0 



(23) 



has been properly defined with the radial variables appearing through the characteristic 
dimensionless ratio 



Mr 3 r i+l 

ieh 



(24) 



As a consequence, the radial path integral, supplemented by the condition r{t) > 0, can be 
given a formal continuum representation 



K l+u (r",r';T)= [ 

J n 



r(T)=r" 



r(0)=r' 



dt 



M 



[r(t)] 2 -V(r(t))\) (25) 



in terms of the usual one-dimensional path-integral pseudomeasure Vr{t). However, our 
subsequent analysis will be based on the time-sliced Eq. (21), which has been called a 
Besselian path integral [16] because of the appearance of modified Bessel functions in the 
pseudomeasure (23). A similar approach can be applied to the energy Green's function 
Gd(?" , r'; E), which through a partial-wave expansion of the form (19) defines a radial 
counterpart 



1 f°° 

G l+U (r", r>- E) = \ dTJ ET ' h K l+u {r'\ r ■ T) . 
i>"> Jo 



(26) 



It should be noticed that Eqs. (21)-(26) explicitly exhibit a property that may be called 
interdimensional dependence [22]: D and / appear in the combination / + v. 
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It is well known that the path-integral expression (21) for the free-particle radial prop- 
agator can be computed analytically and exactly (even before taking the continuum limit, 
N — > oo), with the familiar result [13] 

' Mr'r"" 



r'; T) = V?7> exp 



LU (r' 2 + r" 2 ) 



n+v 



ihT 



(27) 



_2hT 

this is accomplished through the use of Bessel-function identities and recursion relations [23] . 
In addition, the Fourier transform (26) of K^ u (r" ,r';T) [24] provides a closed analytical 
expression for the free-particle radial energy Green's function 



Gff>", r'; E) = -™ y/W 7 l+l/ («r<) Ai +l/ («r>) , (28) 



in which the energy variable appears in the combination 

(29) 

As usual, in Eq. (28), r< (r>) is the smaller (larger) of r' and r", while I p (x) and i^ p (x) are 
the modified Bessel functions of the first kind and second kind (and order p) respectively 

Explicit evaluation of Eq. (21) beyond the free-particle case has been accomplished only 
in a few special cases, including the radial harmonic oscillator [16], using miscellaneous 
limiting procedures in the limit N — > oo. For the case of more general potentials, we will 
next describe a technique that permits the summation of perturbation-theory contributions 
to all orders. 



III. INFINITE SUMMATION OF PERTURBATION THEORY 

In this section we introduce a general technique that consists of a perturbation expan- 
sion followed by its summation to all orders — this is essentially the "perturbation approach" 
pioneered in Refs. [15, 25]. The objective of this combined procedure is to extract "nonper- 
turbative" results, including the emergence of bound states (whenever appropriate), from a 
manifestly perturbative scheme. This is only possible by summing the perturbative series 
to all orders. We generalize the results of Refs. [15, 25] in a number of ways: (i) we present 
the theory in any number of dimensions, for the full-fledged propagator and for the Green's 
function; (ii) we develop it within the time-sliced formulation of path integrals, with empha- 
sis upon possible subtleties involved in nonlinear coordinate changes; and (iii) we explicitly 
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show that the expansion for radial path integrals does not require any additional modifi- 
cations, provided that the Besselian pseudomeasure (23) be used at every order. At the 
very least, this technique will play a crucial role in providing a complete and satisfactory 
derivation of the renormalized physics of the inverse square potential and delta-function 
interactions. 

A. Perturbation Theory 

This procedure is based on the standard resolution of the action in the form 

S [r(t)\ (r", r'; t", f) = S<® [r(t)\ (r", r'; t", t') + S [r(t)\ (r", r'; t", t') 

= S<°> [v(t)} (r", r'; t", t') - f dt V(v(t),t) , (30) 

Jt' 

where [r(t)] (r", r'; t" , t') is the action for a problem whose propagator is already known 
(for example, the free-particle action) while S [r(t)] (r", r'; t", t') provides the additional in- 
teractions V(r, t) in S [r(t)] (r", r'; t", f). 

Even though this beginning step (30) looks like standard perturbation theory, it is in fact 
intended as an exact rearrangement based on infinite summations performed to all orders. 
This resolution can be applied to any path-integral expression: to Eq. (6), for the full-fledged 
D-dimensional problem; and to Eq. (21), for the radial propagator. Just as before, the time- 
lattice versions guarantee that all subtleties are properly accounted for in this formulation. 
Thus, we will start with a lattice action described by the generic Eqs. (5)-(8), as defined by 
the associated partition (3), i.e., 

[ri, • • • r^r-i] (r", r'; t", t') = [n, . . . r^J (r", r'; t", if) + S™ [r x , . . . r N ^] (r", r'; t", t>) 

= S<fW [ri> . . . rjv _j (r - r ' ; t ", - c E V ^*i) ( 31 ) 

3=0 

in lieu of Eq. (30), where S^°' N ^ stands for the corresponding lattice unperturbed action. In 
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each one of these path integrals, the exponential e lS(N) ^ n is to be expanded to all orders, 



exp^SW [r!,...^] (r",r>;t",t') 



'AT-l N-l 

... E-E 

n=0 b'n=° ii=° 



E 



n! 



ih 



ih 



E' 



JV-1 

E 



1 + 



N 



(32) 



"=0 Ul<J2<-<jn=0 

where corrections of order O (1/N) are associated with the terms with repeated indices. In 
order to simplify the notation, while at the same time preventing ambiguities, let us define 



*(o) = — t' t(n+l) = tN — t" 

t(a) = t ja r (a) = r ja (a = 1, . . . , n) . 



(33) 



Then, at each order of perturbation theory, Eq. (32) defines a time ordering with which we 
can associate the partition 

V (n+1) {t',t"} : {fy) =i',*(i), * (2),-.. ,*(„),*(„+!) = *"} (34) 

leading to a resolution of the unperturbed action 

S<M [n, . . . r^] (r", r'; t", t') = [n, . . . r^J (r", r'; t", t') 

n 

= ^ 1 S ( ° ;JV) [ r l>--- r JV-l]( r (a+l)> r (a); i (a+l)^(a)) • (35) 

In Eq. (35) the whole interval T — t" — t' has been divided into n + 1 subintervals T( a ) = 
t( a +i) — t(a), by introducing the points with a = 0, ... ,n + 1, at each order n in 
perturbation theory. Correspondingly, the D-dimensional propagator becomes 



K D (r"y;t",t') = f^K { -\r",r';t",t'), 



(36) 



n=0 



in which the n-th order contribution is 

r r(t")=r" 



K^\v",v';t",t') = lim 

N— »oo 



x e" 



pC)r(t) 

r(t')=r' 
AT-l 

E 

J1<J2< — <jn=0 



1 + 



N 



X II exp<j -5 (0;Ar) [^,...^-1] (r (/3+1) ,r (/3) ;t (/3+1) ,t (/3) ) 

0=0 



(37) 
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with the time-sliced measure displayed in Eqs. (7) and (8). Finally, before taking the limit 
N — > oo in Eq. (37), it is suitable to assume that N 3> n, followed by a rearrangement of 
the time lattice according to the following scheme: 

(i) The whole time interval T = t" — t' is divided into n + 1 subintervals just as in Eq. (35). 

(ii) Each subinterval [i( Q ), i( Q +i)] is partitioned in the usual way into N a parts, just as in 
Eq. (3), i.e., 

^ ^ {^(a) ) t(a+l) } • {^a,jo = t(a) i ta,ji i ■ ■ ■ i ^a,jjv a -i ) ^a,jN a = t(a+l) } • (38) 

(iii) Therefore, the whole interval consists of the following sub-partitions 

pW {t\ t"} = {pW {f , t,} ; {t 1: t 2 } ; ■ ■ ■ ; {t n , t"}} , (39) 

where the notational redefinition implied in Eq. (38) makes the parentheses surrounding the 
subscript a redundant. Thus, the omission of these parentheses amounts to the replacements 



Ha) 



(40) 
(41) 



which will be assumed in our subsequent discussion. With the partition defined by Eqs. (38) 
and (39), a partial measure T>( Na ^r(t) is defined just as in Eq. (7), with the result 



r(t")=r" 



r(f)=r' 



•P v 'r(/) = J} 



a=l 



d D r a 



n 

(3=0 



r(t/3+l)=r/3+i 



r(t/3)=r^ 



(42) 



where the first set of integration factors represents the contribution from the end points of 
each subinterval — conveniently singled out in order to avoid overcounting. Remarkably, in 
Eq. (42), the factors C(M,e,H, D) [from Eqs. (7) and (8)] naturally redistribute in such a 
way that each subinterval inherits the proper measure. 
Finally, in the limit N — > oo, 



N-l 

E 

Jl<j2< — <jn=0 



V(r jn ,t jn ) V(r n ,t n ) 



ih 



ih 



t" rt n 

dt n / dt n -i 
v Jt> 



Jt' 



ih ih 

Then, as a result of the integration measure (42), Eq. (37) takes the form 

KP(r",r';t",f) = f dt n t dt n ^ . . . t dh 
Jt' Jt' Jt' 

V(r a ,t a ) 



(43) 



x 



n 

a=l 



d D r r 



ih 



n^V/m^W^)] , (44) 
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where the contributions from [r(t)\ (rg+i, rp; tp +1: tp) have been converted into the prop- 
agators Kp\rp + i,rp;Tp), as iV — > oo. For the particular case of time-independent Hamil- 
tonians, Eq. (44) can be recast into the alternative expressions 

,d„ V(r a 



kZ\t",t';T) = J] fd D v a 

a=l 

x [ dt n f 

Jt> Jt> 



ih 



dt n -! • • • / dti 



/ dh n [tfgWi, rp; tp +1 - tp)] (45) 
Jt ' p=o 1 J 



and 



/#>(rV;T) = n [d°r a ^ U f°° dTp (rp +1 ,rp;Tp) 5 (t -J> 7 ; 

(46) 

correspondingly, from Eq. (46), the Fourier transform defining the energy Green's function 
yields the result 



G D (r",r';E) = J2G^(r",r';E), (47) 

?1=0 

with an n-th order contribution 

n r poo "| " 

Gg ) (r",r';E) = n / d D r a V(r a ) [ G d\^ + i, TfiE)] . (48) 

a=l L^O J /3=Q 

Direct comparison of Eqs. (45) and (48) shows that the former displays characteristic ex- 
tra integrations with respect to time; these arise when separating the end points of each 
subinterval in the time-sliced expressions. Moreover, the perturbation expansions based on 
Eqs. (36)-(48) admit the usual diagrammatic representations. 



B. Central Potentials 

For the case of central potentials, the same procedure can be directly applied to the 
radial path integrals, with their respective Besselian pseudomeasures (23) at each order. 
Alternatively, Eqs. (36)-(48) could be resolved in hyperspherical coordinates term by term- 
as in the derivation leading to Eq. (19). Either derivation yields the exact infinite summation 

oo 

K l+u (r",r';T) = J2Kl:l(r",r';T), (49) 

n=0 
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for the propagator, where 



n r poo 

K^yy-T) = JJ dr a 



V{r a ) 

ih 



x 



I dt n I dt n -i ... d^Yl r P ] tp+i - tp) 

J t' J t' J t 1 j3=0 



(50) 



Correspondingly, the energy Green's function associated with Eqs. (49) and (50) is 



G l+u (r",r';E) = J2Gtl(r",r';E), 



(51) 



n=0 



where 



a=l 



dr a V(r a ) 



f[[Gil(rp +1 ,rp-E)] . (52) 



0=0 



Finally, in subsequent applications it will prove convenient to introduce the reduced functions 
V(r) and Q^ u (r" ,r'; E), implicitly defined from 

h 2 



and 



G\H(r",r';E) = ~ |V(r")V(r')f^ G\Z{r",r>;E) , 



V(r) 



2M 



2M 



AV(r) 



(53) 



//\-» ,/ /\ 1—1/2 pin) ( ll I. 



(54) 



as well as the sign symbol e ra = sgn (r a ). These substitutions lead to the reduced perturba- 
tion expansion 



g l+u (r",r' ] E) = ^;Ur",r';E), 



(55) 



n=0 



with 



n r pod "| n 

Gi:l(r",r';E) = \"H / dr a e Ta n[C(^ + i,^;£) 

a=l -I fl=o 



(56) 



Notice that, even for n = in Eq. (54), the "unperturbed" function Gi+„(r",r';E) already 
includes information about the perturbation V(r). 



14 



C. Integral Representations 



The perturbation-theory framework of the previous subsection can be further extended 
when V(r) is a monotonic function. This is precisely the case for the class of power-law 
interactions — including the inverse square potential. When this condition is satisfied, the 
sign symbols t Ta in Eq. (56) can be omitted altogether. As a result, the perturbation 
expansion of Eqs. (55) and (56) reduces to the evaluation of 



Q^yy-E) = \ n dr n ... dn 
Jo Jo 

x glH(r", r n , E)gl%(r n , r n _ l5 E) . . . ^(r 2 , r l5 E)g\%{r u r'; E) . (57) 
The formal structure of Eq. (57) can be interpreted to be of the operator form 



#:l(E) = A n [GZiE)\ • (58) 

according to the prescriptions presented in the Appendix, with operator products evaluated 
in real space. This observation leads to a reconstruction of the operator Q^ U (E) in terms of 
the eigenvalues of Q^ U (E), as the following argument shows. 

First, suppose that the eigenvalues and eigenvectors of the self-adjoint operator gf^ v (E) 
are labeled with an index z and known to be i]i +u (z; k) and ^! z> i +v (k) respectively, where n 
is given by Eq. (29). In symbolic notation, 

gl%(E) *,,,+„(«) = rn+viz; k) * z , 1+ „(k) , (59) 

which stands for the integral relation 

dr' Gl%(r, r'- E)V z , l+u (r'; k) = Vl+U (z; K)V z , l+u (r; k) . (60) 

Second, in operator symbolic notation, Eq. (58) is equivalent to the eigenvalue equation 

G\t(E) *,,,+„(«) = A" [ Vl+U (z; K)} n+1 m zHv {n) . (61) 

Third, the corresponding eigenvalue equation for the radial Green's function operator 
gi +v (E) can be obtained by geometrically summing the entire perturbation series, 

g l+u {E) = 1 i - • (62) 

h+Az; k)\ - A 

This equation defines implicitly the operator Q l+V {E). 
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:(o) 



n+l 



f 

Jo 



Finally, Eq. (62) can be recast into the explicit form 



r 

g l+u (r",r';E)= / 



oo 



1 



dz 



h+u(z; re)] 1 - A 



* z ,i +v (r";K)V Zil+v (r';K) , 



(63) 



which stands for the Green's spectral representation (see Appendix). It should be noticed 
that the energy dependence of Qi +V {r" , r'; E) arises from possibly both the eigenfunctions 
ty Zj i +1/ (r; re) and eigenvalues rfi +v (z; re). 

The procedure outlined in the reconstruction of the previous paragraphs illustrates how 
the formal infinite summation can be implemented. However, a caveat is in order: Eq. (63) 
can only be turned into an explicit solution if the eigenvalue problem of Eq. (59) is solved 
first. As an alternative, the solution to Eq. (59) could be directly identified if a corresponding 
spectral representation 



were found. In short, the existence of an integral representation (64) implies directly the 
form (63) for the exact Green's function Qi +V {r" , r'; E). 

For the remainder of this paper we will explore particular applications of this remarkable 
summation technique. 

IV. QUANTUM-MECHANICAL PROPAGATOR FOR THE INVERSE SQUARE 
POTENTIAL 

In this section we start a path-integral analysis of the inverse square potential 



which we conveniently rewrite in such a way that V = 1/r 2 [cf. Eq. (53)] and A be dimen- 
sionless; moreover, A < describes a repulsive potential, while A > describes an attractive 
interaction. Direct evaluation of Eq. (21) leads to complicated integrals and recursion re- 
lations that cannot be easily tamed (unlike the cases of a free particle and a harmonic 
oscillator). Instead, a judicious application of the method of infinite summation of pertur- 
bation theory leads to an exact analytical solution, which will constitute the main result of 
this section. 




(64) 



V{r) 



2M r 2 ' 



(65) 
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Before considering the complete solution, one may notice that Eqs. (27) and (28) describe 
physics of a free particle, for each angular momentum channel, and thereby suggest the fol- 
lowing straightforward heuristic argument. First, one may consider the asymptotic behavior 
of the Bessel function, 

r- _ z . . (z->oo) (p 2 ~ 1/4) 

V2nze Ip{z) ~ 1 , (66) 

in the limit iV — > 00, which according to Eq. (24) amounts to z — > 00. Second, the asymp- 
totics (66) can be applied to each factor in the pseudomeasure (23). Third, each exponential 
factor can be combined with and compared against a corresponding factor arising from the 
inverse square interaction. In conclusion, this shows that the only effect of adding an inverse 
square potential (65) to the propagator (21) is a shift in the angular momentum quantum 
number, 

I + v -> Sl = y/(l + v) 2 - \ . (67) 

Even though this replacement is straightforward in the Schrodinger formulation of the theory, 
its justification in the path-integral approach is far from obvious. Not surprisingly, the 
asymptotic argument has been called into question [6]. Thus, we will first refine the argument 
above and show how to confirm this conjecture (67) in a more satisfactory way. 



A. Integral Representation of the Green's Function for the Inverse Square Poten- 
tial 

The derivation of the result of Eq. (67) can be completed in a number of equivalent ways. 
The most elegant route is based on direct transformations at the level of the energy Green's 
functions. For the inverse square potential (65), the perturbation expansion of Eqs. (55) 
and (56) reduces to the evaluation of Eqs. (57)-(63). 

Remarkably, an integral representation of the form (64) is known, as first pointed in 
Ref. [15]; explicitly, it is the Kontorovich-Lebedev representation [26] 

G[°l(r", r>; E) = -™ ^[W> ' J x+v {kt\ «x'; 1) , (68) 

where 

Jp(Z,m<x) = -5 / dz 1 i K iz (£)K iz (ri) • 69 

7T Jo [Z + V) 
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In Eq. (68) the generalized function J p (£, 77; a) reduces to the product of Bessel functions in 
Eq. (28), 



J p (^rj; 1) = Vwgj® v (ry;E) = /, +I/ («r<)iir, +I/ («r > ) , 



(70) 



when a — 1. 

In the Kontorovich-Lebedev representation [26], which we review in the Appendix, the 



Dirac-normalized complete set of functions ^ z (r; k) = ^J2z sinh(7r,2) / (ir 2 r) K iz (Kr) permits 
the identification of Eq. (69) [supplemented with Eq. (68)] with Eq. (63). Thus, Eq. (59) is 
realized with 



so that the eigenvalues (62) of Qi +U (E) are 

1 



z 2 + (/ + v) 2 ' 



1 



(71) 



(72) 

fa+.C*)]- 1 - A z 2 + si 

It should be noticed that r)i +u (z) is independent of the energy, while ty z is independent 
of (I + v); this is characteristic of the Kontorovich-Lebedev representation. Finally, the 
remarkable reconstruction (63) can be explicitly carried out, 

2M 



G l+u (r",r';E\\) 



h 2 

2 2M 



zsmh(-Kz) 



n=0 

oo 

dz ■ 



7T » J 

2M 



z 2 + s 2 



2 h 2 

V?r"J Sl (Kr",Kr';l) . 



K iz {Kr")K lz {nr') 



h 2 



(73) 
(74) 



Equation (74) is identical in form with (68); thus, the replacement (67) is now proved within 
the path-integral formulation. In short, going back to Eq. (28), 



G l+u (r",r';E\\) = G^ (r",r' ; E) 



2M 
J? 



Vr'r" 7 Si (Kr<)if Si (Kr>) , 



(75) 



while its inverse Fourier transform is identical to K^\r",r';T), so that [cf. (27)] 



K l+U {r"y-T\\) 



M 
ihT 



exp 



%M 
2hT 



( 12 , II2\ 

(r + r J 



Mr'r r 
ihT 



(76) 
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B. Nature of the Solution 



There is only one apparent limitation in the above derivation: the geometric series in- 
volved in the infinite perturbation expansion is guaranteed to converge only for |A| < X\*\ 
where X^ — (I + v) 2 . This condition has been regarded as an actual limitation [15]; how- 
ever, as a restriction, it can be immediately lifted by noticing that the final expression (75) 
provides the desired analytic continuation in the complex A plane. 

Moreover, for each angular momentum state, the radius of convergence of the above series 
is set by the existence of a critical point A = X[*\ at which the nature of the path integral for 
real A changes. In fact, the integral representation (73) explicitly displays that there exist 
two fundamentally distinct regimes with respect to the coupling parameter A (dependence 
which is encoded in si): 

• Weak-coupling regime: A < X\*\ including repulsive potentials. 

The order si of the Bessel functions is real and Eq. (75) has no singularities for real 
k, showing that the system cannot sustain bound states — a manifestation of the scale 
invariance of the potential r~ 2 . 

• Strong-coupling regime: A > X\*\ 

The Bessel functions acquire an imaginary order si = iOi, where 

z = [A-A«] 1/2 , (77) 

and the negative-energy states form a continuous spectrum not bounded from below. 

The pathology displayed in the strong-coupling spectrum can be avoided by the use of 
field-theory regularization and renormalization techniques, as shown in Refs. [9, 10, 11]. 
Accordingly, introducing a short- distance cutoff a, the resulting radial Green's function 
Gi +V (r",r'; E\X; a) inherits the regular boundary condition 

G l+ „(r",r';E\X;a)\ r „ =aoir , =a = 0. (78) 

However, the implementation of Eq. (78) is not straightforward in the path-integral for- 
malism. The main difficulty encountered in a path-integral real-space renormalization with 
the Dirichlet boundary condition (78) lies in the proper implementation of the sum over 
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all paths. A technique for dealing with this difficulty, which was originally introduced in 
Ref. [27], amounts to adding an infinite-strength repulsive delta-function potential cr<5(r — a). 
Therefore, we define the regularized radial energy Green's function Gi +U (r", r'; E\X; a) as the 
limit a — > oo of the more general Green's function Gi +V (r" , r'; E\X; a; a) in the presence of 
the delta-function perturbation. We now turn our attention to this problem. 

V. INVERSE SQUARE POTENTIAL IN THE PRESENCE OF A DELTA- 
FUNCTION INTERACTION 

A. Computation of the Green's Function 

We will now compute the general Green's function G l+l/ (r" ,r'; E\X; a; a) in the presence 
of an interaction 

V^(r) = ~^ + aS{r-a). (79) 

This problem can be most effectively analyzed by treating it exactly to all orders as a 
delta-function perturbation to an inverse square potential. More precisely, the unperturbed 
action includes both the kinetic term and the inverse-square interaction, with their 
exact physics described by Gi+ V (r", r'; E\X), Eq. (75), or alternatively by Ki +U (r", r'; T\X), 
Eq. (76). In addition, the delta-function term V(r) = a5(r — a) is to be regarded as a 
perturbation to which the theory of Sees. II and III applies. 

Then, going back to Eq. (52), each factor J °° dr a V(r a ) is merely reduced to a and carries 
the additional instruction that the replacement r a — > a be made; then, 

G\ll(r",r';E\X;a;<j) = <7 n [G t+V (a, a; E|A)] n_1 G l+lf (r", a; E\X) G t+V (a, r'\ E\X) (80) 

for n > 1, while the term of order zero distinctly remains Gi +U (r", r'\ E\X). The exact infinite 
summation of this series for finite a leads to the familiar result [27] 

Gi +U [r ,r ,E\ X,a,a) = Gi +U [r ,r ,E\X) — — — — — . (81) 

Gi +L ,(a, a; E\A) — L/a 

It should be noticed that Eq. (81) describes the complete physics of a delta-function pertur- 
bation to any known problem with action and exactly described by Gi +ll (r",r'; E\X). 
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B. Applications 



The combined interaction (79) is of current interest, above and beyond its application 
to the regularization of the inverse square potential. Recent work in M-theory has led to a 
particular realization of this combined potential [28] and the use of Green's functions was 
advanced for further analysis [29], within a more general investigation of gravitation in our 
familiar four-dimensional world as arising from a higher-dimensional theory. Even though 
this example corresponds to a very specific model, it shows the usefulness of the examination 
of singular potentials with the techniques introduced in this paper. 

Let us now briefly review the relevant example, which arises in the description of ef- 
fective four-dimensional gravitational effects observed on a (3 + l)-dimensional subspace, 
a 3-brane, embedded in a spacetime with five noncompact dimensions (AdS 5 ) [28]. The 
linearized tensor fluctuations of the metric are subject to a Kaluza-Klein reduction and 
assumed to have a dependence h(x,y) = ip(y)e p ' x in terms of the coordinates x on the 
brane and the extra-dimensional coordinate y. The problem is then reduced to an effective 
one-dimensional Schrodinger equation with respect to the "wave function" i]){y) in terms 
of the extra-dimensional coordinate. Even though direct reduction leads to a complicated 
potential with respect to y, the transformation of variables z = sgn (y) x (e k ^ — l)/k, 
i/j(z) = ip(y)e k ^ y ^ 2 simplifies the effective Schrodinger problem to (h — 1, M — 1) 



--d 2 + 

2 * 



+ <t5(z) 



4>(z) =m 2 ij(z) , (82) 



2{\z\ + af 
with the following specific parameters: 

a = —3k 

a = ± , (83) 

as dictated by the physics of the Kaluza-Klein reduction. More precisely, the additional 
coordinate change £ = \z\ + a, shows that Eq. (82) indeed describes an effective interaction 
of the form (79) with respect to the coordinate £. It should be noticed that the signs in this 
problem (attractive delta and repulsive inverse square potential) are just the opposite of the 
those needed for the regularization of the single inverse square potential, as discussed in the 
next section. 
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The effective dimensionality D — 1 (y — —1/2) of problem (82) implies that the only 
available channels are I = and I — 1, so that I + v — ^1/2 respectively; from this and 
Eqs. (67) and (83), the corresponding order of the Bessel functions is si = 2. Even though 
the repulsive inverse square potential is incapable of producing bound states by itself, the 
additional presence of a one-dimensional attractive delta-function perturbation yields exactly 
one bound state [28]. This is interpreted as corresponding to a massless four- dimensional 
graviton — a fact that has been used to confirm the claim that the experimentally observed 
four-dimensional gravitational fields can arise through a hypothetical scenario of dimensional 
reduction from extra noncompact dimensions. Additional details on this problem and related 
applications are currently under investigation. 

Parenthetically, a related and simple example of the application of these techniques is 
afforded by the two-dimensional delta-function potential V^(r) = —h 2 gS^(r)/2M, for which 
the propagator is also singular and calls for regularization. For example, using a real-space 
regulator a, one may replace 5^ 2 \r) by 5(r — a)/2na. Then, selecting A = 0, v — 0, and 
a = —h 2 g/2M in Eq. (81) and identifying the pole(s), the regularized equation for the bound 
state in the s channel becomes K (K,a) = 2n/g, whose renormalized counterpart agrees with 
the known answer [30], as we have recently shown [31]. 



VI. REGULARIZATION OF THE INVERSE SQUARE POTENTIAL 

We now turn our attention to the problem that was anticipated in Sec. IV: the real-space 
regularization of an inverse square potential with supercritical coupling. In this procedure, 
a short- distance regulator a is introduced and the boundary condition (78) is enforced as 
the a — > oo limit of the Green's function Gi +U (r" , r'; E\X; a; a) in the presence of the delta- 
function perturbation — the second term in Eq. (79). 

Therefore, Eqs. (75) and (81) provide the desired regularized radial energy Green's func- 
tion 

2M y/r'r" 

Gi +v (r",r';E\X;a) = — — — j—^ [ir Si (/ta)J S! (/tr<) - I^K^K^Kr^] K^k^) , (84) 

where k = ^-2ME/h [cf. Eq. (28)]. Equation (84) is in perfect agreement with the result 
from the operator approach to a Green's function: G(r", r'\ E) = u < (r < )u > (r > )/pW[u < , «>], 
where p = h 2 /2M, while u < (r) and w>(r) are the solutions satisfying the boundary condition 
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at r = a and at infinity [32]. 

The poles of Eq. (84), which are implicitly given by 



K Sl (Ka)=0, (85) 

yield the bound-state sector of the theory. The energy levels can be derived by specializing 
to the case when a is small, in which case the small-argument expansion of the modified 
Bessel function of the second kind becomes [10, 33] 



A «- « ' ~ 0> "/wwi sm t e ' ln © " * e J [1 + {z2)] • (86) 

in which <5e ; is the phase of T(l + iOi). In particular, Eq. (86) has an infinite set of zeros 
when the order of the modified Bessel function is imaginary; these zeros are 

Zn = 2 e^i- n7T)/e ' (87) 

[up to a correction factor 1 + 0(z 2 /Qi)], where n is an integer. Furthermore, with the 
assumption that z„ C 1 and 6; > 0, it follows that (—n) < 0, whence n = 1,2,3,.... 
Parenthetically, z n < 1 only if Q t <C 1, so that 5 0; = —76; + 0(Qf) (with 7 being the 
Euler-Mascheroni constant). This argument shows that the energy levels are given by 

2h 2 e^ ( 2im r \ 
E nri = irr~o~ exp --^ , (88) 



Ma 2 r V Oj 
in which n = n r is the usual radial quantum number. 

The regularization of Eq. (85) provides the foundation for the next step: renormalization. 
This final step may be implemented by demanding the dependence of the coupling with 
respect to the regulator; from Eq. (77) this implies the dependence O; = Qi(a) in the limit 
a — > [34], as in Refs. [9, 10, 11]. For example, when this procedure is applied to the 
regularized ground state, the required relation becomes 

2tt 



.,(0) = 



^ + 21n(^) + 2 7 , (89) 



where fx is an arbitrary renormalization scale with dimensions of inverse length and is 
an arbitrary finite part associated with the coupling, such that 

E <«< = - iw exp ^ ■ (90 » 
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The ground-state wave function is obtained in the limit 0. Ja) ^^—^ + , so that, from 
Eqs. (17) and (18), 




The renormalization procedure described above leads to the emergence of an arbitrary di- 
mensional scale — phenomenon known as dimensional transmutation — and violates the man- 
ifest SO(2,l) invariance of the theory [17, 35, 36, 37]. This is an instance of an established 
quantum anomaly, which manifests itself in the three-dimensional molecular realm for the 
interaction between polar molecules and electrons [17]. 

Finally, the scattering sector of the theory can be analyzed in a similar way, with the 
Bessel functions in Eq. (84) replaced as follows: J s (/tr) — > (—i) s J s (kr) and K s (kt) — > 
ni s+1 Hs 1 \kr)/2, and the S-matrix derived with the formulation of Ref. [38]. These results, 
as well as additional properties of the inverse square potential and other singular interactions, 
will be presented elsewhere. 



VII. QUANTUM FIELD THEORY 

Some comments on the quantum field theory case are now in order. The connection to 
this work is twofold: 

(i) The emphasis on singular potentials that we have adopted in this paper can be traced 
fundamentally to their field-theory origin. This explicit connection is illustrated below, in 
the final paragraph of this section. 

(ii) The approach initiated in this paper provides the underlying philosophy and hopefully 
a partial set of ingredients for the more general program of bound states in quantum field 
theory. This generalization is currently in progress. 

In particular, the definitions of Sec. II A can be adjusted as follows. Let us consider the 
path integral 

^0(x,t 2 )=0 2 (x) f ■ \ 

(0 2 (x),t 2 | ^(x),^) = / V<P exp -S[(j>] (92) 

70(x,ti)=0i(x) l n ) 

of the field 0(x) between the states | 0i(x, ti)) and | 2 (x, £ 2 )), corresponding to the prop- 
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agator 



(0 2 (x),t 2 | 0i(x),ii) 



Texp 



*2 



(x)> , 



(93) 



with T being again the time-ordering operator and H the Hamiltonian. For time- 
independent Hamiltonians the "energy propagator" is defined to be 



(M*)\M*))e = 4 (0 2 (x),*2|0i(x),Oe lE(t2_tl)/n 



^0 2 (x) ^-i^ + ^e)" 1 0x(x)^ 



(94) 



(with ie being a small positive imaginary part), which admits a spectral representation 



( fe (x) I ft(x) ) B = E + ^ MM (95) 



E — E n + ie J E - E a + ie 

it 

and permits, in principle, the identification of bound states and scattering states, as well as 
the corresponding wave functionals. 

A subtlety arises in the usual description of bound states at low energies, which requires a 
nonrelativistic (nr) reduction to be performed in one of the following ways: (i) application of 
the nr limit to the original path integral (92), followed by evaluation of the Fourier transform, 
to be subsequently expanded as in (95); (ii) application of the nr limit at the level of Eq. (95). 
It is not obvious that procedures (i) and (ii) are equivalent. In nuclear physics at low energies 
one may wish to integrate out the light field (e.g., pions) and describe the dynamics in terms 
of the heavy field only (i.e., nucleons). Issues about the order in which to take the nr limit are 
also present in this case. Regardless of the route taken, the appearance of singular potentials 
is a generic feature of these nr reductions from field theory to quantum mechanics [39, 40]. 
The proper treatment of these potentials is notoriously difficult and requires regularization 
and renormalization of the corresponding propagator. This particular problem has been 
addressed for the inverse square potential in Sec. VI. 

An intriguing application of these ideas can be found in the recent literature [41] in 
relation to the study of bound states for the Yukawa model, 



(96) 



where \I/ represents fermions and (p the scalar interaction. For two-fermion wave functions 
of total angular momentum J = + , and within a covariantized light-cone technique, the 
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stability of bound states can be assessed by the asymptotics of an effective wave equation 
that simulates a nonrelativistic inverse square potential. Additional subtleties of this prob- 
lem arising from the renormalized strong-coupling sector of the theory are currently under 
investigation. 

VIII. CONCLUSIONS 

In summary, we have formulated the general outline of a program for the evaluation of 
quantum-mechanical propagators and energy Green's functions. We have also thoroughly 
examined the inverse square potential within such path-integral framework. Our analysis 
included renormalization d la field theory in the (strong) supercritical regime, as well as the 
inclusion of an additional delta-function interaction. The general technique as well as the 
particular examples suggest that: 

(i) The treatment of singular potentials — and of bound states in general — is most effi- 
ciently accomplished by the use of the Fourier-transformed or energy Green's function G(E). 

(ii) Infinite summations and resummations in the spirit of perturbation theory capture 
the required nonperturbative behavior associated with the bound-state sector of the theory. 

(iii) Proper analytic continuations are needed in certain regimes. 

(iv) The effective-field-theory program, which leads to singular potentials, requires renor- 
malization in a quantum-mechanical setting — renormalization that, in principle, could be 
implemented with techniques similar to the ones presented in this paper. 

Extensions of this generic program will be presented elsewhere. 
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APPENDIX A: FORMAL OPERATOR STRUCTURE OF THE PERTURBA- 
TION EXPANSION 



In this appendix we summarize the formal algebraic structure of the operator expansions 
presented in Sec. Ill C. In addition, we review the necessary ingredients of the Kontorovich- 
Lebedev transform, properly adjusted to the notation and goals of this paper, and with 
emphasis on the results of Sees. Ill C and IV A. 

A normalized continuous basis for the space of square integrable functions over the in- 
terval X G R is provided by the set of functions {^ z (q)} z£j that satisfy the orthonormality 
condition 

J dq** z „(q)* z ,(q) = 5(z" -z) (Al) 

and the completeness relation 

/ dz^ z (q)%(q') = 5(q-q'), (A2) 
J J 

with the continuous variable z defined over an interval J G R. A particular set {\P ' z (q)} z£ j 
can be constructed by selecting the eigenbasis associated with an appropriate self-adjoint 
operator, with properly chosen boundary conditions. The mathematical theory required for 
such bases, which is extremely subtle, should be studied carefully for every specific case, 
and no attempt is made to develop a rigorous presentation in this paper. Within that 
operational framework, we adopt the usual Dirac notation, which reproduces the required 
formal algebraic structure of Eqs. (Al) and (A2) from 

= 6(J'- J) (A3) 

and 

/ dz \* z ) (*,| = 1 . (A4) 
J J 

In particular, the basis {\q)} provides direct contact with the functions $ z (q) = (q\ty z ). 
These are the functions used in the expansions of Sec. Ill C, where q = r represents the 
radial variable. 

The most important results of Sees. Ill C and IV A rely on the properties of opera- 
tors. With that purpose in mind, an operator A is represented by the "continuous ma- 
trix" A(q",q') = (q"\A\q') in the basis {\q)}. More generally, A admits the representation 
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(\I/ Z //|A|\I/ Z /) in the {\^z)} ze j, in terms of which the sesquilinear expansion 

(q"\A\q>) = [ dz" [ dz'* z „(q") (*AM^z>) ( A 5) 
J J J J 

directly follows by double insertion of the completeness identity (A4). 

Equation (A5) is further simplified for the particular case in which \^f z ') represents the 
eigenbasis of a self-adjoint operator A, with eigenvalues a(z). In that case, 

(q"\A\q')= [ dza(z)* z (q")%(q') (A6) 
J J 

is the required spectral representation of the operator, which we have directly applied, for 
example, to Eqs. (63), (64), and (73). 

The integral transform of the function f(q) by the kernel K(z, q) is defined by 

T {f(q);K(z, q)} (z) = f = J dq K(z, q)f(q) , (A7) 

with appropriate conditions required for the convergence of the integrals involved. Of partic- 
ular interest for the derivations of Sees. Ill C and IV A is the class of integral transforms that 
arise from generalized bases; specifically, given the basis {\^z)} z€ j, the associated transform 
is defined by 

T{f(q);K(z,q)}(z)= [ dqV* z (q) f(q) , (A8) 

J J 

which amounts to K{z, q) = ^* z (q)- Equation (A8 ) is just a particular set of "components" 
of the vector |/), as can be seen with the Dirac notation 

T{f(q);K(z,q)}(z) = (*,\f) , (A9) 

by insertion of the completeness relation for the \q) basis. The abbreviated notation 
f(z) = (z\f) highlights that this class of integral transforms amounts to a mere change 
of basis in Hilbert space — the analogue of the "passive" view of a coordinate transforma- 
tion. Reciprocally, the integral transform (A8) can be used to reproduce the original function 
by means of an inversion identity 

/(g) = / dz*,(q)f(z), (A10) 
J J 

which follows by insertion of the completeness relation for the generalized basis (A2). A few 
of the most common examples of this kind of integral transform are provided by the following 
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partial list of transformation kernels or bases (in ID): (i) exponential Fourier transforms, 
^ z (q) = {2Tx)~ x l 2 e ±tqz \ (ii) trigonometric Fourier transforms, ^ z (q) = a/2/7t r(qz), with 
r = sin or r = cos; (iii) Hankel (Fourier-Bessel) transforms, ^ z (q) = (qz) 1 / 2 J^qz); etc. 
However, this general formal structure is of interest to gain familiarity with more exotic 
cases, such as the Kontorovich-Lebedev transform, which is central to our analysis of the 
inverse square potential in Sec .IV A. 

The Kontorovich-Lebedev representation [26], is based on the modified Bessel function of 
the second kind with imaginary order K iz (q), with the the relevant interval being the real 
half-line X = [0, oo). The notation q = kt proves convenient for direct comparison with the 
equations of Sec. IV A; it is noteworthy that, due to a simple scaling argument (equivalent 
to scale invariance of the inverse square potential), the corresponding Dirac- normalized 
complete set of functions is 



T / \ /22:sinh(7rz) r ^ 

V z (r; k) = \ K iz (kt) . (All) 

V n r 

An important property of the Bessel functions K iz (q) for real variables q and z is that 
Kizio) = K iz (q), i.e., it is real. The orthonormality relation (Al) of the functions 
{^f z (r)} z( zj , with J = [0, oo), amounts to the integral identity 

r°° dr it 2 

/ - K iz „(Kr)K iz ,(Kr) = - , - 8(z" - z') , (A12) 
which can be verified from the integral of the expression [42] 

POO 

K iz n{x)K iz ,{x) = 2 / dt K i{z , +z//) (2x cosh t) cos [(z" - z')t] , (A13) 
Jo 

with respect to x, followed by an exchange of the order of integration and explicit integra- 
tion of the remaining expression [24], with an appropriate analytic continuation. Finally, 
completeness is just equivalent to the Kontorovich-Lebedev inversion identity, 
2 f°° f°° dv 

— dzzsmh(nz)K iz (x)dz -^K %z {y)f{y) = f{x) (A14) 

71 Jo Jo y 

(see Ref. [26]) for functions f(x) defined on the positive real half-line with appropriate 
conditions required for convergence. 
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